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1. Introduction 

Green functions (GFs for short) of colour singlet local QCD currents have proved to be an 
essential tool for a successful description of hadronic physics. Confinement binds quarks 
and gluons into colour singlet hadrons and hence it makes little sense to compute matrix 
elements with those fundamental particles as asymptotic states. The process of hadro- 
nisation takes place at an energy scale of the order of Aqcd and thus it is essentially 
non-perturbative. QCD currents can be written as 

Yj l/^j^^^y- = Yg(a;):, (1-1) 

where F and A" are spin and flavour matrices, respectively, and A denotes a (summed) 
colour index. The F matrix stands for 21, 2^75, 7^^, 7'^75 and a^"^ in the case of scalar, 
pseudoscalar, vector, axial-vector and tensor sources, respectively.^ The currents Jp(x) 
qualify as interpolating fields for mesonic states in the sense of having non-vanishing 
matrix elements between the vacuum and the hadronic state: ( Hj | Jp | ) 7^ 0. This 

"""In our conventions 75 — 170717273 and a'"' = ^ [7'', 7"]- 
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permits to obtain matrix elements of hadronic states from GFs through the LSZ reduction 
formula. Let us define the two- and three-point GFs in momentum space: 

Uf,{p) = ild'xe^P-^0\T{jnx)4m\0), 

^f2UP,q) = d4xd4ye*(P-+''-^)(0|T{Jf(x)j|(y)J|(0)}|0). (1.2) 

These definitions can be trivially generalised for more currents. 

QCD exhibits a spontaneous breakdown of the chiral symmetry due to quantum effects. 
It is believed that the operator responsible for this phenomenon is the quark condensate 
(qq), which acquires a non-vanishing expectation value. Phenomenologically, this feature is 
reflected through the appearance in the spectrum of the so called pseudo-Goldstone bosons. 
At very low energies (long distances) E <^ nip, GFs can be computed as a perturbative 
series in the (small) momentum carried by the current and the quark masses in the so 
called Chiral Perturbation Theory (xPT) which constitutes a low-energy effective 

field theory (EFT). In this formalism the non-perturbative effects are encoded in the a 
priori unspecified low-energy constants (LECs) that grow in number as we go further in 
the chiral expansion. 

At high energies (short distances) the QCD currents in the GF approach each other 
and one can expand their product as a tower of local operators in the Wilson operator 
product expansion (OPE) In momentum space this is tantamount to an expansion in 
inverse powers of momenta, only being valid in the deep Euclidean region —p^^ m^. The 
vacuum expectation values of these local operators are known as vacuum condensates and 
encode the non-perturbative effects js], which start to become relevant at intermediate 
scales. Then we have an explicit separation of the short distance effects in the Wilson 
coefficients and the long distance effects in the condensates. The arbitrariness of this 
splitting is reflected in the dependence of both the Wilson coefficients and the condensates 
on the renormalisation scale and the employed renormalisation scheme. 

At intermediate energies the most reliable tool is an expansion of QCD in powers of 
[0) H) 1]) being the number of colours of the QCD gauge group SU(A''c). This 
region is populated by hadronic resonances driving the strong dynamics. At leading order 
in 1/Nc there must be an infinite number of such states for any given set of quantum 
numbers. GFs of a finite subset of these resonances can be obtained as tree-level diagrams 
of an effective Lagrangian with resonances as explicit degrees of freedom, known as RxT 



|1C, 11|. Nevertheless, we lack a method to handle an infinite tower of such states, and 
then the so called minimal hadronic ansatz (MHA) is usually used, meaning that we only 
consider the lowest lying resonances. One is then able to perform a matching of the three 
regimes and thus to estimate the values of LECs present in the xPT Lagrangian [|l^, 

The power of EFTs is enhanced by matching onto the more fundamental theory in a 
region where both descriptions are sensible, and running down with the renormalisation 

^Throughout this paper we will assume the chiral limit of QCD, that is we will set the masses of the 
light quarks u, d and s to zero. 

^An approach in the framework of Pade approximants to approximate GFs in the large- A'^c limit which 
encompasses saturation with a finite immber of resonances, can be found in Ref. |14]. 
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group equation (RGE) to a scale where the EFT is useful. This results in a resummation 
of large logarithms that would otherwise spoil the perturbative expansion. As hadronic 
theories, the matching between xPT and RxT is straightforward since up to a certain 
mass scale the spectra of both EFTs coincide, and the matching can be accomplished 
be formally integrating out the heavy degrees of freedom. However, the situation changes 
drastically when matching RxT onto QCD, because the spectrum no longer has any particle 
in common. The reason for it is that we are going through the chiral phase transition. Then 
one is forced to compare GFs computed in the different EFTs in momentum regions where 
their domains of validity overlap. 

The matching relation is better understood for GFs which are order parameters of 
chiral symmetry breaking (order parameters for short). These GFs have a vanishing Wilson 
coefficient for the identity operator (that is, the purely perturbative contribution) to all 
orders in as in the chiral limit, and hence its leading contribution stems from the {qq) 
condensate. Since this operator is also responsible for the chiral symmetry breaking those 
GFs encode essential information on its mechanism. 

This method has been used for two-point GFs including several orders of Og corrections 
|15], but only at leading order for three-point GFs [16, |l^, 19, ^ (see also |21] for 



a matching with two multiplets of vector-meson resonances). It is the purpose of this 
paper to push the matching up to 0{as). Unfortunately, the 0{as) corrections to the (qq) 
Wilson coefficient C^gg^ for three-point GFs are not known, and we will concentrate on 
their computation, relegating the details of the matching to a forthcoming paper. Finally, 
for the case of the Hvt GF, we will show that we can perform the matching to the one-loop 
OPE result with a single multiplet of vector-meson resonances. 

The paper is organised as follows: In Section Q we discuss the renormalisation scale 
dependence of the GFs and how this affects the matching; in Section ^ we compute the 
0{as) corrections to the OPE two-point GFs and match the (VT) result onto R^T with 
only one multiplet of vector-meson resonances; in Section Q we outline the calculation of 
the 0{as) corrections to the OPE three-point GFs and argue what the difficulties are when 
matching onto the hadronic representation; in Section |5| we present our conclusions. 



2. Scale dependent Green functions 

A strong motivation for computing the 0{as) corrections to C^gg) is to become sensitive 
to the renormalisation dependence of the quark condensate in full QCD. As is well known, 
since a QCD current involves two quark fields in the same space-time point, it is not 
enough to renormalise the quark fields to get a finite result. The current itself must be 
renormalised in addition: 

J^ = ZrJr, (2.1) 

which defines the renormalisation constant Z-p- The superscript B denotes /i-independent 
bare currents while currents without superscripts correspond to renormalised (hence gener- 
ally /i-dependent) ones, where fi signifies the arbitrary renormalisation scale. It is custom- 
ary to express the //-dependence through the anomalous dimension 7r, which is defined 
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as 

_ fi dZr dJr , , 

7r = ^ -7— , - fJ- -7— = ir Jr ■ (2.2) 

Zr a// d/i 

The anomalous dimension depends on the couphng = Os/vr, and in perturbation theory 
has an expansion 

7r {as) = 7?^ + al + ^^^Kl + . . . . (2.3) 

At leading order in one can easily calculate the anomalous dimension of a current 
with the following master formula (of course, anomalous dimensions are gauge invariant): 



= ^r-^7'^7-r7.7M, (2-4) 

where the algebra should be performed in four space-time dimensions. Plugging in the 
structures for the different F's one finds that vector Vff and axial- vector A"; currents have 
zero anomalous dimension. This is a general result that stems from the fact that in the 
chiral limit both currents are conserved. In the MS scheme which we follow in this paper 
the renormalisation procedure is mass independent and so the result for the anomalous 
dimensions holds for finite quark masses as well. Similarly, one can show that the quantities 
niq-S'^ and niq-P"' are (j. independent, if the currents are normal-ordered. Whence it follows 
that 

IS = IP = —-r^ = -Im, 2.5 
niq d/i 

to all orders in ctg. The /i dependence of the renormalised current is reflected in the GFs 
themselves. We can write the dependence through a RGE ^ 



o n 

a \ - 



ni...„ = 0. (2.6) 



Here, 7^ are the anomalous dimensions corresponding to the currents appearing in ni...„. 
QCD vacuum condensates in general also have a non-vanishing anomalous dimension, 
and in particular for the quark condensate {qq), as should be clear from the above, 7(gg) 
coincides with 75 of Eq. (2.5).^ 

Let us now consider the case of the 3-point Green function 11555 with three scalar 
currents. In the OPE it gets its first contribution from {qq), and at leading order its 
Wilson coefficient is //-independent. On the RxT side the situation is a bit different. It 
is well known that each scalar and pseudoscalar current insertion in the chiral theory is 
accompanied by a {qq) factor, so [19| 



liTss = C§f{p, q) {qq)if^) , H^g = U{p, q) [{qqKf^)]' , (2.7) 



''This argument is spoiled if the GF needs subtractions, that is it does not correspond to a physical 
quantity. Since we are deahng with order parameters, our GFs do not need any subtraction. 

^There exists a sophistication in that for a proper separation of long and short distances in the OPE, 
non-normal-ordered condensates should be employed, but then m{qq) is not RG-invariant p^ . 
However, here we are working in the chiral limit, and for vanishing quark masses this problem is absent. 
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and on first sight, we cannot match one onto the other because the fi dependences appear 
different. What is happening is that RxT includes ah ag orders in a non-perturbative 
fashion, while in QCD the coefficient function C^-^^ must exhibit some /i dependence to 
account for the missing scale factors. This situation easily generalises to other Green 
functions, and we can write a RGE for the coefficient function C^qq)(/x), such that Eq. 
is fulfilled: 



d 



Now the scale dependence of the Green functions in both theories of Eq. 
in principle the matching could be performed. 



(2.8) 
agrees and 



3. Two— point GFs and matching 

To prepare the discussion of the 3-point functions below, we begin with reviewing the 
calculation of the ag corrections for C^gg) in the case of two-point order parameter GFs. 
In this simpler scenario, we will discuss the appearance of infrared (IR) divergences at 
intermediate stages of the calculation and the renormalisation of {qq). For 2-point functions 
in which the quark condensate is the leading contribution in the chiral limit, there are only 
two GFs to be considered: H^p and Ily!^". 

The fact that vector and axial-vector currents are conserved in the chiral limit has 
implications on the GFs, known as Ward identities. For the case of the vector-tensor GF 
they imply that Ily^{p) = 0, and so taking into account the antisymmetry of the tensor 
current, together with parity, we can parametrise it as 

i^vT)t,,ip) = i6'''{g^.uPp-g,pPu)IlvT{p'). (3.1) 

For the H^p GF, Ward identities have somewhat deeper consequences and in fact com- 
pletely fix it in the chiral limit: 

{UAp)f{p) = 2i6^'^{qq), (3.2) 

where a detailed derivation can be found in Appendix A. Again, we see that the scale 
dependencies on both sides of Eq. ( |3.2| ) are identical since the anomalous dimensions of the 
scalar current 75 and of the quark condensate agree. Furthermore, one finds that in the 
chiral limit Hap is saturated by one pion exchange. 

Let us now concentrate on the Hy^" GF. The relevant diagrams are shown in Fig. |l|. 
Diagrams (c), (e) and (f) are infrared divergent since they involve gluons attached to quark 
lines with zero momentum. The same diagrams contribute to Hap and Eq. ( p.2| ) shows 
that it is free from IR divergences. So they cancel when adding the six diagrams, and the 
same occurs for the Ilyj' GF. Since these divergences cancel at the end we might choose 
any method to regularise them. In what follows, we will adopt dimensional regularisation, 
which simplifies computations notoriously. In this scheme diagrams (e) and (f) are zero: 
they are scaleless and convert the IR divergence of diagram (c) into an ultraviolet (UV) 
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one. In particular, for this diagram we obtain 



n 



ic) 
VT 



Cf {qq) J , N 
— ^1 -(3 + a) 



(3.3) 



where a denotes the gauge parameter.^ As we see, by itself this diagram is still gauge 
dependent, which means that the other diagrams are required in order to obtain a gauge- 
invariant result. 

Summing up all diagrams, we find a remaining divergence of the form — Cp ag/e. Part 
of this divergence (namely — 3C_pas/4e) is absorbed in the renormalisation of the quark 
condensate, that is {qq)^ = Zi^gg^{qq){fj,), while the remaining divergence is removed with 
the renormalisation of the tensor current, = ZtJt{i^), leading to 



1 



1 



(3.4) 



From the renormalisation constant Zt^ we can also derive the first coefficient of the tensor- 
current anomalous dimension, 7^^^ = Cpjl^ in agreement with the simple formula of 
Eq (p.4|). For the full vector-tensor Green function at the next-to-leading order, we then 
find 



OPE 

VT 



{qq) 



1 + ttsCp 



log 



p2 



1 



2Nc 



which is of course gauge invariant. For Nq flavours Cp 
mation corresponds to the large-iVc limit. With our result of Eq. 
the RGE for the coefficient function is trivially satisfied: 



+ Oialp- 



(3.5) 



^f- where the approxi- 



|3.5|), we can check that 



d 

/" t;- + 7m + 7T 



C. 



VT 



d 

/i — +2asCF 



VT 

m 



(3.6) 



®A11 calculations in this paper have been performed in an arbitrary covariant gauge. The dependence 
on the gauge parameter a cancels in all our final, physical results, constituting a good check for them. 
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Let us now match this result onto the large-A'^c prediction of R^T. Since the {VT) 
GF in the chiral hmit does not need to be subtracted, it is completely determined by its 
spectral function. Then it can be regarded as an observable and thus we can directly match 
its expression in different approximations, such as OPE and RxT. However, we still have 
the problem that the tensor current requires renormalisation, and thus, the (VT) GF is 
scale and scheme dependent. This renormalisation dependence would then be reflected in a 
scale dependent coupling of the tensor current to vector mesons fyifJ.) on the hadronic side. 
Since we prefer to work with hadronic quantities which are explicitly scale independent, 
another possibility is to multiply /y (/i) by an appropriate scale factor Rt{h), which results 
in a scale independent tensor decay constant fy. This is analogous to the definition of scale- 
invariant 5-factors, which parametrise hadronic matrix elements of four-quark operators, 



in the case of weak hadronic decays |25]. Therefore, we define 



fvif^)RT{fJ') = fvip) exp 



(2) 

Tt 
(3i 



IT {as) 



da. 



as{fj.) + 0{a 



(3.7) 



Nf 



' 27 



|Ia.(^) + 0(a2) 



The anomalous dimension of the tensor current is known up to order |26, 27 1, and thus 
one could even extend Eq. (^^)- However, at the order considered here this does not make 
sense, since we only stay at the next-to-leading order level. Now multiplying our result 
( |3.5| ) for the Green function with the scale factor i?j'(/i), it is a trivial exercise to convince 
oneself that Rt{i^) ny^^(p^, /x) is scale independent at the considered order. Nevertheless, 
it should be kept in mind that it still depends on the renormalisation scheme, for example 
on the scheme in which (qq) is renormalised. 

In principle, since the next-to-leading order result for the (VT) GF contains a log- 
arithm in the dynamical variable p^, let us strongly emphasise that an infinite tower of 
resonances would be required for a sound matching to the RxT. Still, as a simple-minded 
approach, we will next consider the aforementioned minimal hadronic ansatz (MHA). This 
amounts to the assumption that a single resonance is enough to correctly describe the 
physics in a certain energy regime."^ Tensor sources in chiral Lagrangians were first intro- 

m 



duced in Refs. |29|. In Ref. |Q the matching for the (VT) GF was performed at 0{a^), 
and here we will include the 0{as) corrections. Again multiplying the GF with the scale 
factor in order to obtain a scale-invariant quantity, the hadronic ansatz reads: 

fv fv '^V 



(3.8) 



The precise definitions for the decay constants fv and fy can be found in Eq. (53) of 
Ref. ^ 



Eq. (p^) is in principle assumed to be valid at all energies at leading order in 



^For a discussion of the subtleties arising when matching the fuU tower of vector-meson resonances to 
the OPE for the (VT) GF the reader is referred to Ref. S. 
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1/Nc, since it incorporates chiral symmetry and the correct high-energy behaviour. It 
can be expanded in inverse powers of p^, permitting a direct comparison with the OPE 
in Eq. ( |3.5| ). However while Eq. ( |3.8| ) is expUcitly scale independent, ( |3.5D contains a 
logarithm which compensate the running of the tensor source and the quark condensate. 

To perform the matching in practice, we choose a particular matching point and scale. 
First of all, to sum up the logarithm, we will employ the scale fj? = —p^ = M^. Then, 

should be large enough so that only keeping the first term in the OPE is a good 
approximation, while it should not be too large so that only putting one resonance on the 
hadronic side is reasonable. From these considerations, we would conclude, that M should 
be in the range 1-2 GeV. For the matching relation, we then find 



fv fv '^v 



■2[a,(M2)]-i^{A}(gg)(M2) 



1 



985 
486 



8357 
11616 



Nc } UsiM^ 



(3.9) 



where in the curly brackets, we have also included the numbers corresponding to the large- 
Nc limit. Eq. ( |3.9| ) can be viewed as a refinement over the analogous estimate of Ref. pl| . 

Let us finally come to a numerical analysis of Eq. (^]^) . Employing the central values 
fv = 221 MeV and My = 775 MeV |^l|, as well as the value for the quark condensate 
{qq){2GeV) = - (267 MeV)^ pi, and Qs(Mz) = 0.119, we obtain 



fj 



195 ± 55 MeV , 



(3.10) 



where the quoted uncertainty dominantly results from a variation of the matching scale M 
in the range 1-2 GeV, and to a lesser extent from either taking the renormalisation-group 
coefficients in full QCD, or the large-A'^c' limit. The large matching-scale dependence 
of our result reflects the imperfection of the matching. At a scale of 1 GeV, the scale 
dependent vector-meson tensor coupling reads: /y(lGeV) = 167 ±47 MeV. Given the 
large uncertainties from the matching scale, this finding is in surprising agreement to the 
leading order result (1 GeV) = 165 MeV [^] and to determinations of the tensor coupling 



from QCD sum rules and lattice QCD of Refs. pOl 



4. Three-point GFs 

Let us now concentrate on the main topic of this work: the corrections to C^qg) for 
the three-point order-parameter GFs. Three-point GFs are of great interest because of 
several reasons. First, unlike the two-point case, there is a quite large amount of them 
that are order parameters of the chiral symmetry breaking. Second, in the framework of 
RxT they involve vertices between resonances and so they are useful for studying how the 
resonances interact. Third, there are many 0{p^) LECs that can be determined with 
these GFs And fourth, by means of the LSZ reduction formula we can relate the GFs 
with form factors entering the calculation of many interesting hadronic observables. Some 
phenomenological applications of these Green functions can be found in Ref. |^, ^] and 
references therein. 

There are eight such GFs and we will sort them according to their intrinsic-parity 
and Lorentz tensor rank. For rank-two GFs we can use the chiral Ward identities and 
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the transformation properties under parity for splitting them in terms of only one or two 
Lorentz tensors (made up of g^'^ , e^^"^ ^ and the external momenta). This amounts to a 
great simplification on their calculation. When two or more currents have the same Dirac 
structure one can also exploit Bose symmetry. We already saw that the flavour structure 
of two-point GFs is trivial: it is proportional to a Kronecker delta (there is no other SU(3) 
rank-two tensor). For the present case there are only two rank-three SU(3) tensors, the 
totally antisymmetric structure constant /"^'^ and the totally symmetric d"'^'^ tensor. Either 
one or the other will appear as a global factor depending on the transformation properties 
of the GF under time reversal. 

The diagrams giving rise to the gluonic corrections are shown in Fig. |2|. Diagrams (h), 
(i) and (j) are analogous to those in Fig. |l| (c), (e) and (f), respectively. Diagrams (d) 
and (e) are IR and UV safe, even though they have gluons attached to zero momentum 
quark lines. Much as we did in Section ^, we regularise those divergences in dimensional 
regularisation, renormalise the quark condensate as in the last section and after summing 
up all diagrams IR divergences will cancel out. The rest of the (UV) divergences are 
absorbed in counterterms as usual, resulting in a finite (but scale dependent) result. 

Loop corrections manifest themselves as logarithms, dilogarithms and constant pieces. 
In general we will have the following decomposition: 



{<!<}) 



Cf 



Lp log 



) +Lglog (--^ 



+ Lr loe 



+ LdCo + Lc 



(4.1) 



where r = — {p + q) and the Li are /^-independent meromorphic functions of the squared 
external momenta. This simple structure arises because in the chiral limit all internal lines, 
*In this work we use the convention e°^^^ = +1. 




Figure 2: Gluonic corrections to C(qq) for three point GFs. 
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either quark or gluon, are massless. The massless, scalar three-point integral Co{p , q 
its all dilogarithms and its ex 

1 

7x 



collects all dilogarithms and its explicit expression reads [34|: 

d^fc 1 

(2^ k'^{p+k f{q-k)^ 

^ i Ll2 5 ^— ^ - Ll2 



A-g2 


— p2 


_|_ j,2 


A + g2 


+ r2 




A-g2 




+ p2 


A + r2 




-q' 


A-r2 


— p2 


+ g2 



+ Li2 5 5— — ^ - Li 





— p2 


+ r2 


A + g2 


+ p2 











A + g2 


_|_ j,2 


— p2 


A-r2 


-P' 


+ 9^ 


A + r2 


+ p2 


-g2 



Li2 ( I - L.a j I . (4.2) 

where A is the well-known Kallen function A (p^, q'^,r'^) = (p^ + 9^ — '^^)^ — Ap^q^. Up to 
order a^, the ^ dependence of C^gg) then corresponds to 

d 

- ^ C'(gg) (/^) = as— (^p + i^g + ^r)- (4.3) 



From the general expression Eq. (4.1) we see that there is no possible choice of /i that makes 
all logarithms cancel simultaneously. Then we are forced to include an infinite number of 
resonances in order to match onto RxT. Moreover, besides logarithms we also have the (//- 
independent) dilogarithms, and it is unclear how to match those to the sum of resonance 
contributions. Let us next explore the different sectors. 

4.1 Zero rank GFs: (SSS) and (SPP) 

Before starting the discussion of this sector we remind the reader that GFs of the type 
[SSP) or [PPP) are forbidden by parity invariance of the strong interactions. The Lorentz 
structure of this sector is trivial, as it has no Lorentz index; time reversal symmetry fully 
determines the flavour structure: 

n5S5(5PP)(/>9',r2) = d'^''''Iisss{SPP){p\q\r^) , (4.4) 

Bose symmetry requires Ilsss{p'^,q'^,^'^) to be totally symmetric in its three arguments 
and also Uspp{p'^, g^ r2) = nspp{p'^,r^, q^). 

The total anomalous dimension of these GFs is 7 = 875 = — 87^, as discussed in Sec- 
tion 1^. Thus when calculating them in RxT (or xPT) one finds that they are proportional 
to {qq)^-^ On the other hand in the OPE side the first non-vanishing contribution comes 
from the (single) quark condensate. The RGE for its Wilson coefficient, Eq. ( |2.8| ) is then 

d ^ \ ^sss{SPP), 



^^-27™JCJ:P^^^^(m) =0, (4.5) 



®In xPT with its hadronic degrees of freedom, there is no dynamical generation of the QCD renormal- 
isation scale fi {i.e. it does not come from the UV divergences of Feynman diagrams) and thus it is fully 
contained in _Bo = — {qq)/F^ , the quark masses and unphysical couplings like Hj and //J. Furthermore, all 
LEG of the additional operators in the chiral Lagrangian which arise in the presence of tensor sources 
also depend on fi. The QCD renormalisation scale jj, is not to be confused with the chiral renormalisation 
scale fj,^, showing up when computing chiral logarithms. 
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where the Wilson coefficients C^^^^^^^\fi) have been defined as 



^SSSiSPP)ip',q',r^fi) = Cj^f(^^^)(p2,g2^r2,/x)(g-<z)(/.) . (4.6) 



At 0{a^) their expressions read: 



OfS = - ^ Mf^ '!^ •■^) . = ^ [P* - ie - r^l ^ (4.7) 

Our exphcit results for the 0{as) corrections to these expressions (and for the remaining 
Wilson coefficients) are rather lengthy and thus have been relegated to Appendix B. 

In Ref. 1 19 1, the leading order expressions ( |4.7| ) have been matched to hadronic repre- 



sentations which took into account constraints arising from xPT and RxT, thereby allowing 
to estimate some of the LECs in these frameworks. Going to the next-to-leading order, 
now we have to face two complications. First of all, we have the question of scale depen- 
dence of the GFs. This problem could be treated like the {VT) GF of the last section: 
by defining quantities which are explicitly scale independent through multiplication with 
appropriate scale factors. (This problem cannot even be addressed at the leading order.) 
Secondly, due to the presence of logarithms (and dilogarithms) the matching definitely 
requires an infinite tower of resonances, making it much more involved. For this reason, we 
postpone a detailed discussion of the matching procedure for three-point order-parameter 
GFs to a forthcoming publication.^'^ 

4.2 Odd intrinsic parity sector: (VVP), (AAP) and (VAS) GFs 

These GFs are rank-two Lorentz tensors and so in principle they can be written in terms of 
many independent Lorentz tensors built from external momenta, the metric tensor g'^'^ and 
the Levi-Civita symbol ei^'^'^l^. However chiral Ward identities imply ^123 = ^123 ~ ^• 
Besides, the transformation properties of the odd-intrinsic-parity sector can only accom- 
modate a Levi-Civita symbol contracted with the two external momenta. Together with 
time reversal invariance we are left with the following structure: 

{llvVPiAAP)[VAS]T';;(.P,q) = I^VVPiAAP)lVAS]{p',q\r')e^uaf3p''q''d'^'''(.d^'l[r'"']. (4.8) 

Bose symmetry implies ^vvp{AAP){p^t<1^t''''^) = ^vvP{AAP)iQ'^ ^p'^ t^^)- The total anoma- 
lous dimension of these currents is 7 = 7s = — 7^, and this means that both in Rx? (xPT) 
and leading order OPE the result is proportional to a single quark condensate. Moreover, 
the RGE for the Wilson coefficient Eq. ( p.^ ) reduces to the fact that it is /^-independent: 

d ^VVPiAAP)lVAS] _ „ . , 

Analogously to the (SSS) and (SPP) GFs, we then define the Wilson coefficients as 

nvyp^AAPWAs]{P^q'y,^.) = c;^7(^^^)[^^^](p^g^r2) (mm . (4.10) 



^"An additional complication in satisfying quark-counting rules that requires the matching of higher- 
point functions with an infinite set of resonances is discussed in Ref. [p5|. 
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At lowest order they are given by 



C, 



VVP 

m 



+ + 
q2 p2 J.2 



c- 



AAP 

m 



p2 _j_ g2 



q2 p2 J.2 



c, 



VAS 



p- _ 



j2 p2 J.2 



(4.11) 



The next-to-leading order corrections are again presented in Appendix B. Here a tech- 
nical comment is in order. Since we work in dimensional regularisation, there is a question 
how to treat 75. In all our computations we have employed a fully anticommuting 75. 
Either two 75's appear in a trace and can be cancelled before taking the trace, or, in the 
odd-parity sector, we can first perform the 7 contractions before taking the trace, and are 
then left with traces of only four 7-matrices and a 75 which are unambiguous. As expected, 
in the odd-parity sector our Wilson coefficients turn out to be scale independent. 

4.3 Even intrinsic parity sector: (VVS), (AAS) and (VAP) GFs 

These last GFs are also rank-two tensors and so they can be built out of many Lorentz 
structures as well. The chiral Ward identities are not so simple as in the odd-intrinsic- 
parity sector and read: 

\abc , 



p''iqniiivvsr;:{p,q) 

p^{qn{nAAs)f:{p,q) 



0, 

2 {qq) d'^'^' % 



p2 



{IIVAP)Z iP^ <l) 



q" {UvapTJI^ (p, q) 



2{qq) r''^ 



(4.12) 



where we have also used time reversal invariance. These results rely on the fact that the 
{AP) GF is fully determined by the chiral Ward identity Eq. (|3.2|) , and they must be 
satisfied in any sensible description of the strong interactions. For instance, concentrating 
on {AAP) and (IMS'), in RxT this means that once we contract with one momenta, only 
the pion pole can survive. In the OPE they also have deep implications: since there is 
no other condensate than {qq) in Eq. ( 4.12| ), it implies that the contribution of higher 
dimension condensates vanishes when contracting with one external momenta, and since 
there is no Ug factor all contributions to the {qq) beyond leading order must also vanish 
when the contraction is performed. Eq. ( |4.12[ ) allows us to write the GFs as 



(n 



VVS) 



V abc 



2 2 2\ 
J ,p ,r ) 



abc 



P'^'ip, q) J'vvsip^q^r^) + Q^'ip, q) Qvvs {p\ 



KAsT'^'il'yy) = P'''iP,q)^AAs{p',q'y)+Q'''ip,q)GAAs{p',q'y) 



+ 2{qq)^ 
q'^ p^ 



P^%P, q) TvAP {p\ q\r^) + Q^^P, q) QvAP {p\ g^ r^) 
{p + 2q)^q^ g^"" 



-'i{qq) 



^2 ^ 2 



(4.13) 
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where the transverse tensors P^'^ and Q^'^ are defined as 



Q^''{p,q) ^ p^q'^q' + q^p^p" -{p-q)p^q'' -p^^g'^" . (4.14) 

T and Q are two scalar functions satisfying Bose symmetry in the first two arguments for 
the (yVS) and {AAS) GFs: 

^[Q]vvsiAAS) {P^q^r^) = ^[G]vvsiAAS) {q',P^r^) ■ (4-15) 

At lowest order the determination of the scalar functions J- and G is straightforward, 
but once we go beyond this level their direct computation turns out to be rather com- 
plicated. Instead it will be much simpler to concentrate on the determination of linear 
combinations of these factors, obtained by taking appropriate contractions in Eq. ( 4.13| ) 



(i.e. those contractions that do not reduce to a Ward identity of Eq. ( [4.12] )). They take 
the form 

A 



KVS = ^{p' + q'- r') Tyys - + ?>p\')QvvS 



qi^Pu^vvs = - ^ ^vvs + ^ (p^ + - ?'^) Gvvs , 

9^^u ii^As = ^ (p^ + 9^ - r^) ^AAs - + ^p'^q^)GAAs - {qq) ^^2q2 ^ ' 

-nUV '^'T- ^ / 2 2 2\^ /— \ i.P^ ~^ 

Qf^P'^^AAS = -4-^AA5+g(p +q - r ) Qaas + {qq) , 

9.^ Kap = 1{P' + ^'- r') TyAP - + ^p\')QvaP + m ^^'\'fr''^ 



q^p, = --JVAP^\ + q'- r') QvAP - (qq) ^^2^2 '^^ " (4-16) 

Once these contractions are known, we are in a position to determine the scalar functions 
and G by inverting Eqs. ( |4.16| ). The contractions have the same symmetry properties under 



exchange of momenta as J- and G, due to Bose symmetry. Eqs. ( 4.16 ) admit an expansion 
in as, and the terms proportional to (qq) contribute only at zeroth order in as- Thus, 
for the radiative corrections we are concerned about, the third and fifth lines reduce to 
the first, while the fourth and sixth reduce to the second. The same discussion concerning 
RxT and the OPE as in the odd-intrinsic-parity sector applies also here, and again Wilson 
coefficients turn out to be scale independent. We define those Wilson coefficients as 

9^ll'^WvS{AAS)[VAP]y^^) - '-'{qq) yP '9 jW^/lMjj 

^.P'^KvsiAAswAP^i^^) = cl^r'^^'^^'^^%\q\r^)m{,). (4.17) 
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The leading-order results for these coefficients are found to be: 



vvs 



1 



(99) ^l rpl 

^gAAS _ 1 



(2 2^2,/2| 2\ 2 r)4 

[p -q ) + [p +q )r -2- 



r 



ill) pi cp- r2 



qpAAS _ ^J^/x I »2 2\ 



Like before, the complete results of our computation at order ag can be found in Ap- 
pendix B. 

5. Conclusions 

In this work we have made a step forward in our understanding of quark-hadron duality, 
by computing order ag corrections to the leading quark condensate contribution for 2- and 
3-point Green functions which are order parameters of the chiral symmetry breaking. We 
have concentrated on these order-parameter GFs, since in this case the matching between 
the OPE result and the corresponding GFs obtained in a low-energy effective theory like 
xPT or RxT is more transparent. 

The matching further simplifies in the large-A'^c? limit of QCD, because then, the 
hadronic spectrum only consists of an infinite set of zero-width resonances. Still, generally, 
at the next-to-leading order, the appearance of logarithms in the kinematical variables 
requires the inclusion of the full infinite number of those resonances in order to be able to 
reproduce the logarithms present on the OPE side. 

For the 2-point {VT) GF, we have discussed the matching in more detail. The inclusion 
of as corrections allows to trace the dependence of both, the OPE and the hadronic side, 
on the short-distance renormalisation scale ^i. These dependencies show up in the quark 
condensate as well as possible anomalous dimensions of the initial currents on the OPE 
side, while in the low-energy theory, besides the condensate, the LECs of the tensor sources 
explicitly depend on ^ [|2^. To be able to do some numerics, finally we have assumed the 
minimal-hadronic-ansatz as a crude approximation, and have matched the {VT) GF with 
a single multiplet of vector-meson resonances. This allowed to estimate the vector-meson 
tensor coupling fy. 

The main and most involved aim of our work was the calculation of the Ug corrections 
for the leading OPE term of 3-point order-parameter GFs. Details of the computation for 
the different types of 3-point GFs to be considered have been discussed in Section 4, while 
the rather lengthy final results have been relegated to Appendix B. 

Additional complications for the matching in the case of 3-point GFs arise from two 
facts: firstly, due to the presence of three independent kinematical variables, no choice of 
the renormalisation scale allows to resum all of the logarithms which are present in our 
results. Secondly, in addition in the next-to-leading order result also dilogarithms arise 
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and it is unclear how to match them even with an infinite set of resonances on the hadronic 
side. Therefore, we postpone a detailed discussion of the matching in the case of 3-point 
GFs to forthcoming work. 
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Appendix A: non— renormalisation of the (AP) GF 



In this appendix we demonstrate that in the chiral limit Ward identities completely fix the 
{AP) GF, to be fully saturated by single pion exchange. Lorentz and SU(3) invariance 
imply 

{UAp)f{p) = i6'''UAp{p')p,. (A.l) 

From the definition of time-ordered product and the conservation of the axial-vector cur- 
rent it follows that 



d'^ r{^^(x) P\0)} = 6ix") [AS(x), P%0)] . (A.2) 
On the other hand Fourier transformation relates derivatives with external momenta 

o) = 0. (A.3) 



J d^x e^-f (9^ + ip^) ( I T{^;i(x) P^{0)} 



The right hand side of Eq. (A.2) can be easily computed since the Dirac delta function 
forces the time dependence on both fields to be the same: 



<5(x°-yO) \A-,ix),P\y) 



— 6^'' S{x) + d^^" S'ix) } (x - y) , (A.4) 

Uf J 

where S{x) = Ylnf 'If If flavour singlet scalar current. Then from Eqs. ( A.l ) to (|A.5D 

it follows that 

2i 



ah J2 ■ 



5''^(0|5(0)|0) = 2i6'"'{qq), 



ab 



Uf 



(A.5) 



where we have used the fact that vacuum is flavour blind. We can identify Hap = 2 {qq) /p^ 
and then demonstrate Eq. ( |3.2D . 



Appendix B: explicit expression for the Wilson coefficients 

In this appendix we present the complete expression for the corrections to the quark 
condensate Wilson coefficients. We have split our results as shown in Eq. (f4.l|), in terms 
of coefficients multiplying logarithms, dilogarithms and polynomial terms. The results can 
be found in Table |^. 
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{SSS) 



Lp 



4 



[ 4/ + + - 6 g2 r2 - (g2 _^ r^)] 



4 ^4^4 + g4 ^^4 _ 6g2^2 _ 3^2 (^2 ^^2^j 



P 9 



[p^ + 4 



6 1 »,6 



2p4 (-^2 _^ ^2-) _ 2 g4 (-^,2 _^ ^2-) _ 2 j,4 ^^2 _^ ^2^ j 



(- - - 5r4 + 14p2 r2 + 14g2p2 + 14^: 



Lc 

Lp 

L, 
Lf 

Ld 

Lr 



[r^ + — Ap^ — (r^ + q"^) — Qr"^ q^ 



4 [4g4^^4_^4^3^2(^2_^2)j 



_^[4r4 + g4_/ + 3r2(p2_g2)] 

[ ^6 + ^6 _ ^6 _ 2p2 (^4 + ^4) + 2 - g2 ,2) (^2 + ^2)] 

(- 5r4 - 5g4 + + I4r2g2) 



-2 
-2 



p6 _|_ ^2 p4 _|_ p2 ^2 ^2 _j_ q,2^ _ 2 (|^2 _j_ ^,2^ ||^2 _ ^,2^ 2 
^6 _j_ ^2 g,4 _j_ p2 q2 ^2 _|_ ^2^ _ 2 (r^ + p^) (r^ — p^) ^ 



-4 

_^p2g2j,2 



[ 2 (/ _ 3 p2g2 ^ ^4) ^4 _ 2 ^2 (^2 + ^2) (^4 + ^4) 
2r6 (p2 + ^2) + ^8 + (p2 _ ^2)2 ^^4 + ^4 _ ^2^2) 
(p2+g2+4^2^ 



-2 

_^p2q2^2 

-2 
^p2g2^2 



p6 _ 3 ^2 p4 + p2 ^2 (5 ^2 + g2) ^ 2 (r^ - g2) 3' 
^6 _ 3^2 qA_^p2q2 ('5^2 _^p2^ _^ 2 (r^ -p2)3 



\p'2q'2r2 



4 r*^ + 6 (p2 + g2-) _ ^,2 ^^2 _^ 3^2^ ^^2 _^ 3 ^2^ _^ ^^2 _ ^2-) ^ ^^2 _^ ^ 



Table B: corrections to the Wilson coefficients for three-point GFs. 
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-4 


^ [2r6 (p2 + g2^ _^ gp2^2^4 2r2(p2 


+ (p^ + q") 










— r^ + {p^ — q^) ^ [p"^ + — p^) 






Lc 






-^{f+q^ 4r2) 






(VAS) 


Lp 


-2 


"/ + 3p4 (2g2_r2) _p2q2 (5^2+^2) +2 








-2 r 

\p'2q'2y.2 \_ 




+ (g2 + 6r2)p4_ (654 + 6^4_5^2^2^^2_^6_^2r6-g'^r2] 




[p6 + 7^2p2(g2 _p2) _ ^6 _ 4^6 + 2 (Sp^ + ^2) ^4 _ 3 (p4 _ ^4) ^2] 


Ld 


Xp-^q 


4 


(3 + 2 g4 + 6 ^4 ^2 _^ 2 r^) p2 _ 


^-g2)2(^4 + ^4) 








+ p*^ - p6 (3 g2 _^ 2 r^) - 2 r2 g2 p4 j 






Lc 














Lp 






_^[2(g4-.4)_p4+^2 (,2_ 






Lq 






_^[2(p4-r4)-g4 + g2 (,2_3p2)] 










^;r^ [6p2 g2^4^4_^4_^4_^2(^ 


' + q')] 




Ld 






_^[2r6 {r^+p^ + q^ 3^^ (^2 + ^2)] 




Lc 






[3r2(p2+52_^2)_4^2^2] 






Lp 






^ [g2 (2g2_2p2+^2) _3^2 (p2 






Lq 






^ [p2 (2p2_2g2^^2) _3^2 (^2 












_^[3(/ + g4)_3^2(^2 + ^2)_2p2^2] 




Ld 




2 

— T 

p'^q^ 


[p^ _ p4 ^2 _|_ ^^2 _ j.2^2^^2 _|_ ^2-j _ ^2 ^2(^^^2 _|_ j,2-jj 




Lc 






^(p2 + g2-2r2)A 









Table B: Ug corrections to the Wilson coefficients for three-point GFs. 
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Lp 



2[qi -r^f -p^+p^ (5r2-3(72 



2 ip^ - 



[2r6+ (-3r4+/ + g4_3^2g2^ + ^2^] 



[3^4 _ ^2 _ 3^2 ^^2 _^p2^j 



Lr 



p^r^ 



_^[^2 (g2_3p2 + 3,2)+2g2 (g^-p^)] 
^ [r^ (p2 _ 3^2 _^ 3^2^ _ 2p2 (-^2 _p2^j 
_^ [_ 2p2 ^2 + 3 (^4 + ^4) _ 3 ^2(^2 + ^2)] 
_2_^ _ 3^4 ^2 ^ 3p2 ^4 ^ (^2 _ ^2)3j 

1 (^2 + ^2 + 2r2)A 



Lp 

Lq 
Lf 

Ld 

Lr 



_^[2{q'-r^)+p'-p^ {5r^ + q')] 



,2^21 



piqifi 



2 ■ [p2(g2_^4^2^ _2p4_g4_2^4_^g: 
_^[p4_^4 + 4,4 + ,2 (p2_^2)] 

- 2 g2 / + g4 _ ^6 _^ 2 - (3p2 + g2^ j 



Lq 

Lf 

Ld 

Lr 



3 (p2 + 3g2 
^(3/ + g2-.2^ 



.2 ^2^ 



r 



^5 — T 

pZqz 



p^{q^-p^) {q^+p^-r^) 

p6_3^2^4^^2^2 (2^2 + 3^2) _ (^2 _ ^2) 2 + ^2) 



^[p^-q^ + 2r^}X 



Table B: Os corrections to the Wilson coefficients for three-point GFs. 
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